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ABSTRACT 

Intermediate Mass Black Holes (IMBHs; IO^ '^^^Mq) are thought to form as relics 
of Population III stars or from the runaway collapse of stars in young clusters; their 
number and very existence are uncertain. We ran N-body simulations of Galactic 
IMBHs, modelling them as a halo population distributed according to a Navarro, Frenk 
& White (NFW) or a more concentrated Diemand, Madau & Moore (DMM) density 
profile. As IMBHs pass through Galactic molecular/atomic hydrogen regions, they 
accrete gas, thus becoming X-ray sources. We constrain the density of Galactic IMBHs, 
f2,, by comparing the distribution of simulated X-ray sources with the observed one. 
From the null detections of Milky Way Ultra-Luminous X-ray sources, and from a 
comparison of simulations with unidentified sources in the IBIS/ISGRI catalogue we 
find a strong upper limit n, < 10-'^nb{< lO^^rib) for a DMM (NFW) profile, if 
IMBHs accrete via ADAF disks. Slightly stronger constraints (il, < lO^'^fif, for a 
DMM profile; O. < 10"^^^^ foj. ^ NFW profile) can be derived if IMBHs accrete with 
higher efficiency, such as by forming thin accretion disks. Although not very tight, 
such constraints are the most stringent ones derived so far in the literature. 

Key words: black hole physics - methods: A^-body simulations - Galaxy: general - 
X-rays: general 
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1 INTRODUCTION 

There are strong observational evidences for the existence of 
two classes of black holes (BHs): stellar mass BHs, with mass 
from 3 to 20 Mq (Orosz 2003), thought to be the relics of 
massive stars, and super massive black holes (SMBHs) in the 
mass range W^~^Mq, hosted in the nuclei of many galaxies. 
Recently, the existence of a third class of Intermediate Mass 
BHs (IMBHs) has been inferred. They are characterized by 
masses in the range from 2OM0 to a few xW^Mq (see van 
der Marel 2004 for a review). 

Several IMBH formation mechanisms have been pro- 
posed: (i) IMBHs could be the relics of very massive metal 
free stars (Heger & Woosley 2002), (ii) they could form 
in young clusters via runaway collapse of stars (Portegies 
Zwart & McMillan 2002), or (in) could have been built up 
in globular clusters through repeated mergers of stellar mass 
BHs in binaries (Miller & Hamilton 2002). Some recent ob- 
servations indicate that IMBHs could exist in the core of 
globular clusters (Gebhardt, Rich & Ho 2002, 2005; Gerssen 
et al. 2002; van den Bosch et al. 2005). Their number is 
nearly unknown. In principle, IMBHs could contribute to 
all the baryonic dark matter (van der Marel 2004): their 
density Q, could be essentially equal to that of luminous 
baryonic matter, Qb. lum = 0.021 (Persic & Salucci 1992; 



Fukugita, Hogan & Peebles 1998), and equal to 50% of all 
baryons {Qb = 0.044, Spergel et al. 2003) . Only weak con- 
straints on the IMBH mass have been inferred from dynam- 
ical studies of the Milky Way. For example, the observed 
velocity dispersion of stars in the Galactic disk requires that 
halo BHs masses are < 3 x 10^ Mq, if the Milky Way dark 
halo is entirely made of compact objects (Carr & Sakellar- 
iadou 1999; see also Lacey & Ostriker 1985; Wasserman & 
Salpeter 1994; Murali, Arras & Wasserman 2000 and ref- 
erence herein). Other dynamical constraints on the IMBH 
mass can be derived by imposing that they do not disrupt 
too many Galactic globular clusters (Moore 1993; Arras & 
Wasserman 1999). By this request, Klessen & Burkert (1996) 
found that, if the dark halo of the Milky Way is exclusively 
composed by IMBHs, their mass must be < 5 x 10* Mq . For 
the same reason, the halo BHs cannot represent more than 
2.5-5 per cent of the dark halo mass, if they are as massive as 
10^ M© (Murah et al. 2000). However, constraints on IMBHs 
obtained from globular cluster disruption are very uncertain, 
as we do not know what are the characteristics of globular 
clusters when they form, and how many of them have been 
destroyed. It could even be that IMBHs have played a role in 
determining the current number and distribution of globular 
clusters (Ostriker, Binney & Saha 1989). 

In this paper we explore an alternative way to dynam- 



2 M. Mapelli, A. Ferrara & N. Rea 



ically constrain the IMBH density, based on the proposed 
identification of Uitra-Luminous X-ray sources (ULXs) witli 
IMBHs. ULXs are, by definition, point sources with X-ray 
luminosity higher than 10^^ erg s^^, exceeding the isotropic 
Eddington limit for a 10 Mq BH (sec Colbert & MiUer 2005 
for a review). ULXs have not been found, up to now, in 
the Milky Way; but they are present in many spiral and 
starburst galaxies (Swartz et al. 2004). ULXs tend to be as- 
sociated with star forming regions; but they often lie near, 
not in them (Mushotzky 2004). 

ULXs were initially identified with SMBHs with a low 
accretion rate; but this interpretation was found to be in 
conflict with their position in the host galaxies, far off 
from the galaxy center (Colbert & Miller 2005). Later on 
ULXs have been suggested to be high-mass X-ray binaries 
(HMXBs) with beamed X-ray emission (King et al. 2001; 
Grimm, Gilfanov & Sunyaev 2003; King 2003). Even though 
low luminosity ULXs {Lx ^ 5 x 10^^ erg s~^) are consis- 
tent with this HMXB scenario, the highest luminosity ULXs 
show various characteristics which can be hardly reconciled 
with the beaming model (Miller, Fabian & Miller 2004) , such 
as the existence of a ionized nebula surrounding some bright 
ULXs (Pakull & Mirioni 2002; Kaaret, Ward & Zezas 2004). 

An intriguing hypothesis, at least for these highest 
luminosity ULXs (Miller et al. 2004), is their identifica- 
tion with accreting IMBHs. This idea is also supported by 
some observational evidences. First, the spectra of many 
high luminosity ULXs have a soft component well fitted 
by a multicolor black-body disk, whose inner temperature 
is typical of BH masses in the IMBH range (Miller et al. 
2004; Colbert & Miller 2005). In addition, high luminos- 
ity ULXs often show long term variability on timescales of 
months to years (Kaaret et al. 2001; Matsumoto et al. 2001; 
Miyaji, Lehmann & Hasinger 2001) and quasi periodic oscil- 
lations (Strohmayer & Mushotzky 2003), inconsistent with 
the beaming scenario. King & Dohncn (2005) propose that 
very high luminosity ULXs in interacting galaxies can be 
IMBHs hosted in the merging satellite and whose accretion 
is activated by tidal forces. 

Many studies have been dedicated to check the pos- 
sibility that ULXs arc IMBHs accreting in binary sys- 
tems (Baumgardt et al. 2004; Hopman, Portegies Zwart 
fe Alexander 2004; Kalogera et al. 2004; Portegies Zwart, 
Dewi & Maccarone 2004; Hopman & Portegies Zwart 2005; 
Patruno et al. 2005). However, the observed population of 
ULXs is not well reproduced by this binary system sce- 
nario (Blecha et al. 2005; Madhusudhan et al. 2005), mainly 
because the ULX phase of simulated IMBHs is too short. 
A better agreement between simulations and observations 
can be obtained only by considering very massive IMBHs 
( ^ 1000 Mq ; Baumgardt et al. 2005) . In addition, very few 
optical counterparts have been detected for ULXs so far and 
can unambiguously be identified as companion stars (Liu et 
al. 2005; Colbert & Miller 2005). Then, it is stiU open the 
possibility that ULXs are IMBHs accreting gas during the 
transit through a dense molecular cloud, as recently sug- 
gested by Krolik (2004) and by Mil & Totani (2005). 

This paper is aimed at exploring in detail this last hy- 
pothesis by dedicated N-body simulations (Section 2). In 
particular, we simulate a Milky Way model and we derive 
an upper limit of the density of IMBHs, by requiring that 
no ULX is produced in the Milky Way by IMBHs passing 



through molecular clouds (Section 3). Next, we study the 
distribution of both ULX and non-ULX sources produced 
by IMBHs passing through molecular clouds (Section 3) and 
atomic hydrogen regions (Section 4). We finally compare the 
derived distributions with observations (Section 5). 



2 NUMERICAL SIMULATIONS 

The simulations have been carried out using the parallel N- 
body code Gadget-2 (Springel 2005). The simulations were 
performed using 8 nodes of the 128 processor cluster Avo- 
gadro at the Cilea (http:/ /www. cilea.it). Our aim is to gen- 
erate a suitable N-body model of the Milky Way, in which 
we embed a halo population of IMBHs. 



2.1 Milky Way model 

To reproduce the Milky Way wo simulated an exponential 
disk and a Hernquist spherical bulge, whose density profiles 
are given, respectively, by the following relations (Hernquist 
1993): 

Pd{R,z) = exi>-(R/Rd)sech^(z/zo) 
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where Md (Mb) is the disk (bulge) mass, Rd is the disk 
scale radius, zo is the disk scale height, a is the bulge scale 
length and r = \J R? -\- . We choose a = 0.2 Rd, consistent 
with Kent, Dame & Fazio (1991; o = 0.7 ± 0.2 kpc). The 
value of 20 is quite difficult to constrain. Recent observations 
(Larson & Humphreys 2003; Yoachim & Dalcanton 2005) 
suggest the presence of two components in the thin disk of 
the Milky Way: a "young star forming" thin disk (with scale 
height ~ 200 pc) and an "old" thin disk (with scale height 
~ 600 pc). Then, we assume 2o = 0.1 J?d = 350 pc, which 
is approximately an average of the scale height of these two 
components and correlates in a simple way with Rd- 

Disk and bulge are embedded in a rigid dark matter 
halo, whose density profile is (Navarro, Frenk & White 1996, 
hereafter NEW; Moore et al. 1999): 



Ph{r) = 



(r/r.)'' [l + (r/r,)"]('3-^)/-' 



where we choose (a,/3, 7) = (1,3, 1), and pa = 
being the critical density of the Universe and 
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(3) 

Per it Pcrit 

(4) 



3 ln(l-hc)-(c/(l-hc))' 

where c is the concentration parameter and rs is the halo 
scale radius, defined by Ts — R200 / c; -R200 is the radius en- 
compassing a mean ovordensity of 200 with respect to the 
background density of the Universe, i.e. the radius contain- 
ing the virial mass M200. Given the concentration c and the 
Hubble parameter^ H{z), R200, M200 and the circular ve- 
locity at the virial radius, V200, are related by the following 



^ We adopt H{z ^ 1) = Hq = 71 km Mpc"^ in agreement 
with first year WMAP results (Spergel et al. 2003) 
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expressions. 
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Table 1. Initial parameters for the Milky Way model. 



-R200 — 



M200 = 



WH{z) 

WGH{z) 



(5) 



G is the gravitational constant. In our simulations we fix 
c — 12, V200 = 160 km s^^ (sec Klypin, Zhao & Somerville 
2002), yielding M200 = 1.34 x 10^^ M©, 7?200 = 225 kpc and 
rs = 19 kpc (Table 1 reports the initial parameters) ; finally, 
we use for H{z) its actual value Hq. 

Rigid halos can induce m = 1 instabilities in the disk. 
To check this effect, we performed test simulations with 
a non-rigid halo (with halo particles ten times more mas- 
sive than disk particles); we did not observe significant dif- 
ferences in the evolution with respect to fixed-halo mod- 
els. Since simulations with non-rigid halos are prohibitively 
time consuming for the very high resolutions required by 
the problem (see Section 2.3 for details), we have chosen to 
adopt a rigid halo. 

To derive Md, Mt and Rd we followed the prescriptions 
of Mo, Mao & White (1998) , imposing that the disk is a thin, 
dynamically stable and centrifugally supported structure, 
whose mass is a fraction of the halo mass and whose angular 
momentum is a fraction of the halo angular momentum. In 
particular, our best, stable model is obtained for a choice of 
the spin parameter A = 0.035, where A = J \E\'^^^G~'^ M''^^^ 
(J, E and M being the angular momentum, the total energy 
and the mass of the halo, respectively). Requiring that Md + 
Mb « 0.04 M and that Md : M;, = 4 : 1 (in agreement with 
Kent et al. 1991; Freudenreich 1998; Binney & Merrifield 
1998), we obtain M<i = 4 x 10^° Mq and = 1 x lO" Mq. 
Our choices are in agreement both with the best model of 
Milky Way described in Klypin et al. (2002) and with the 
COBE measurements of the bulge mass (Af[, = 1.3 ± 0.5 x 
10^° Mq, Dwek et al. 1995). For these values, the scale radius 
of the disk is Rd = 3.5 kpc, consistent with recent estimates 
(Binney & Merrifield 1998). Given the uncertainty on the 
Md/Mt ratio, we also made some test simulations for A4d '■ 
Mb = 5:1 observing no significant differences in our results. 
In order to account for the SMBH in the nucleus of the Milky 
Way, we located a point mass Msmbh = 3.5 x 10^ M© (Ghez 
et al. 2003; Shodel et al. 2003) at the center of the rigid halo. 

Initial velocities of disk and bulge particles are simu- 
lated using the Gaussian approximation (Hernquist 1993) 
for dispersion velocities. This choice introduces a transient 
behavior, represented by outwards propagating rings of over- 
density from the warmer disk center, as it was already noted 
by Kuijken & Dubinsky 1995 (see also Kazantzidis, Magor- 
rian & Moore 2004; Widrow & Dubinski 2005). In agree- 
ment with the findings of Kuijken & Dubinsky 1995, in the 
highest resolution runs (when the mass of each particle is 
m ^ 10^ M0 and the total number of particles approaches 
one million) this transient is stronger; nevertheless, it always 
disappears within about 1 timescale^, when the system re- 
laxes into a new equilibrium configuration. We consider this 
new relaxed configuration as initial condition for our analy- 
sis. This procedure is legitimate in our case, since we are not 



^ The timescale of our simulation is defined as the rotation period 
of the simulated galaxy, i.e. about 0.27 Gyr. 
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investigating processes such as disk instabilities, but we are 
only interested in the dynamical evolution of halo IMBHs. 
After relaxation, we continue the simulation for about 5 Gyr, 
i.e. from redshift z w 0.5 until today, about half of the time 
elapsed from the last major merger (Governato et al. 2004). 
This allows us to follow the evolution of an already relaxed 
and nearly unperturbed (by mergers) Milky Way. During the 
entire simulation disk and bulge remain perfectly stable. 



2.2 Intermediate mass black holes 

How many IMBHs are hosted in the Milky Way ? What 
is their spatial distribution ? These axe yet unanswered 
questions. Nevertheless, we need an Ansatz on the IMBH 
number, mass and distribution to generate the initial condi- 
tions of our simulations. A reasonable estimate for the ini- 
tial IMBH number follows from Volonteri, Haardt & Madau 
2003. Assuming that the IMBHs are born mi>a fluctuations 
(with 1/ = 3 — 3.5) collapsing at a given redshift, Volonteri 
et al. (2003) derive the density of IMBHs at the formation 
redshift Q,, / as 



fl,,f = [1 - erf (i^/V2)] 



M(v) 



(6) 



where \i — eri{v / \/2)\Q,m is the fraction of the Universe mat- 
ter in halos with M > M(v) {yim = 0.27 being the matter 
density), as derived from the Press & Schechter (1974) for- 
malism, and m,,//M(i/) is the fraction of mass of the halo 
collapsed in IMBHs {rn,j being the average initial IMBH 
mass, and M{v) the mass of the va peak halo). For exam- 
ple, assuming that IMBHs form in 3 cr fluctuations collapsing 
at redshift z = 24, the corresponding halo mass is M(3) = 
1.7 X 10^ Mq (Barkana & Loeb 2001). Under these assump- 
tions equation (6) gives il.j — 10"'* fit (m._//10''MQ). 

Given fi,,/, one can roughly estimate the number of 
IMBHs in the Milky way, N, as 



N. 



fib 



in.,f 



(7) 



where Mb^MW = (0.5—1) x IO^^Mq is the mass in baryons of 
the Milky Way. Adopting a value of / = lO"** Q,b, we find 
A^. ~ 10*. Instead, if we assume that IMBHs form in 3.5 cr 
fluctuations collapsing at z = 24, this number becomes N, = 
5 x 10^. We will adopt equations (6) and (7) to calculate how 
many IMBHs to include in our simulations. 
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How massive are the IMBHs today? We have assumed 
that their average mass at formation was m,j — 10^ Mq. 
However, it is likely that they accreted gas for some period 
of their life (Ricotti & Ostriker 2004; Madau et al. 2004; 
Shapiro 2005; Volonteri & Rees 2005). The duration and 
the efficiency of accretion are highly uncertain, making hard 
to determine the amount of accreted mass. Shapiro (2005) 
suggests that the IMBH mass evolution m, (t) , assuming Ed- 
dington rate accretion, can be written as: 



m,{t) — m,j exp 



1 



t 



iSalp 



(8) 



where e is the radiative efficiency (e ~ 0.1) and tsaip is 
the Salpeter time {tsaip ~ 0.45 Gyr). In our case t is the 
fraction of IMBH life during which it accretes at the Ed- 
dington rate, i.e. t = fdutyhirth, where hirth is the time 
elapsed from the IMBH formation (« 13.5 Gyr) and fduty 
is the fraction of turth during which the IMBH accretes. 
Assuming fduty ~ 0.01 (fduty ~ 0.03 for quasars ed, z ^ 6, 
Steidel et al. 2002), we obtain m,{t) « 10 m,,/. For our 
choice of m,j- = 10^ M©, this means that the average mass 
of IMBHs today is m.(t) « 10** M©, consistent with the 
value 1.8 x 10* M© of the recently detected IMBH candidate 
in the globular cluster Gl (Gebhardt et al. 2005) and with 
previous theoretical estimates (Volonteri et al. 2003). This 
estimate is affected by a number of uncertainties, and we 
consider it only as an educated guess. 

Due to accretion, the current density of IMBHs fi. will be 



ion., / = lO"^f7(,(m.,//lO^A'/0), (9) 



our reference value. 

We note that other models predict very different val- 
ues for Q,. For example, Salvaterra & Ferrara (2003) de- 
rived Q, « 0.1 r^i,, under the assumption that Population 
III stars are the sources of the observed near-infrared excess 
with respect to galaxy counts. One of the aims of this pa- 
per is to check which part of the SI. range is allowed by the 
link between IMBHs and ULXs (see next section). For this 
reason, we also carried out two runs adopting the estimate 
by Salvaterra & Ferrara (2003). 

The last problem we have to address is the selection of 
initial conditions for the position and velocity distribution of 
IMBHs. White & Springel (2000) suggested that remnants 
of Population III stars should be much rather concentrated 
inside present-day halos. N-body cosmological simulations 
by Diemand, Madau & Moore (2005, hereafter DMM) seem 
to support this idea; they also show that the present spatial 
distribution of objects formed in high-cr fluctuations depends 
only on the rarity of the peak in which they are born. In par- 
ticular, DMM flnd that the spatial distribution, in present 
halos, of objects formed in a i^cr fluctuation is well fltted by 
a modified Navarro- Frenk- White profile: 



p.(r-) 



Ps 



{rjr^)^ (1 + (r-/r,)°')(f--^)/° 



(10) 



where p^, a and 7 are the same as defined in the previous 
section; = Vsj fv is the scale radius for objects formed in a 
1/ (7 fluctuation (with = exp(i//2)), and = 3-1-0.26 1'^ ''. 
As DMM simulations are coUisionless, they cannot take into 
account the possible formation of binaries containing IMBHs 
(eventually with the central SMBH) and the occurrence of 
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Figure 1. Density profile of IMBHs for a DMM (case Al; solid 
line) and for a NFW distribution (case A2; dotted line). 



three-body encounters, which likely lead to the ejection of 
one of the involved IMBHs (Volonteri et al. 2003). Thus, the 
actual IMBH distribution could be slightly more "diluted" 
than that obtained by DMM. 



2.3 Description of runs 

We made three different sets of simulations. A, B and C, 
whose characteristics are described in Table 2. Simulations 
labeled as A are characterized by S7. = 10~'^f2t, (correspond- 
ing to IMBHs formed in 3 a fluctuations), simulations B 
have f2. = lO'^f^i, (corresponding to the Salvaterra & Fer- 
rara 2003 model), and the simulation C has SI. = 3 x 
(corresponding to IMBHs formed in 3.5 a fluctuations). Sim- 
ulations Al, Bl and C adopt the DMM spatial distribution, 
assuming that IMBHs form in 3 ct (Al, Bl) or 3.5 a (C) 
peaks. Instead, runs A2 and B2 were performed assuming 
that IMBHs follow a normal NFW profile. In Fig. 1 we com- 
pare the two considered distributions of IMBHs, i.e. DMM 
and NFW. In each simulation we evolved about 10^ parti- 
cles, each having a mass of 5 x 10'* A/0 (included the IMBH 
particles). These particles are divided in ~ 8 x 10^ disk par- 
ticles and ~ 2 X 10^ bulge particles, plus a variable number 
of IMBH particles: 2000 in each run of the group A (cor- 
responding to 10* IMBHs of 10* Mq ; see equations (7) and 
(8)), 2x10^5 in each run of the group B (corresponding to 
lO'^ IMBHs of 10* Mq) and 100 in the run C (corresponding 
to 500 IMBHs of 10* Mq). CPU time limits require that we 
consider only equal mass particles, with mass no lower than 
5 X 10* Mq (included the IMBH particles), making impos- 
sible to investigate the impact of dynamical friction on the 
IMBH spatial distribution. 
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Figure 2. Snapshots at t = 1.4 Gyr (about 5 Galactic timescales) of the simulated Millty Way; only 1/200 of the total bulge and disk 
particles are plotted. Big dots indicate the IMBHs passing through the molecular disk in the case Al. Left panel: particle positions in 
the X — y plane; central: in the x — z plane; right: in the y — z plane. 



Table 2. Initial parameters for the IMBHs. 



Lx{p,v) = ri c'G^ mI pv ^, 



(11) 



Run 


Number of IMBH particles 




IMBH profile 


Al 


2000 


10-3 


DMM° 


A2 


2000 


10-3 


NFW'' 


Bl 


2x10^ 


10-1 


DMM 


B2 


2x10^ 


10-1 


NFW 


C 


100 


3 X 10-5 


DMM 



"Diemand, Madau & Moore 2005. 
''Navarro, Frenk & White 1996. 



3 IMBHS ACCRETING MOLECULAR GAS 

As discussed in the Introduction, one of the possible expla- 
nations for the existence of ULXs is that they are IMBHs, 
accreting both in binary systems (Patruno et al. 2005) or 
in molecular clouds (Mii & Totani 2005). Here we investi- 
gate the possibility that ULXs are IMBHs accreting gas dur- 
ing their transit within a molecular cloud. We also checked 
how many non-ultra-luminous X-ray sources {Lx < 103^ erg 
s"^) could be produced by IMBHs passing through molecu- 
lar clouds. 



3.1 IMBH density constraints from ULXs 

The X-ray luminosity3 of a BH with mass M, , expected from 
the Bondi-Hoyle accretion in a gas cloud, can be expressed 
as (Edgar 2004; Mii & Totani 2005): 



where r) takes into account the radiative efficiency and the 
uncertainties in the accretion rate, c is the hght speed, p 
the density of the molecular cloud and v — (u^ -f al^^j -f- 
Cs)^''^! being the IMBH velocity relative to the gas; omc 
and Ca are the molecular cloud turbulent velocity and gas 
sound speed, respectively. Recently, Krumholz, McKee & 
Klein (2006) have shown that accretion rates in a turbulent 
medium might slightly differ from the above one. Because of 
the many uncertainties in our model, we do not attempt to 
deal with these subtleties. From equation (11) and following 
Agol & Kamionkowski (2002), Mii & Totani (2005) derive 
the number of ULXs with luminosity higher than Lx as**: 
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erg s~ 



(12) 



where /i is the mean molecular weight {/i ~ 1.2-2.3 depend- 
ing on the fraction of H2 molecules), A^, is the number of 
IMBHs in the Milky Way (see equation 7) and 77 is the ra- 
diative efficiency. The correct value of rj is completely uncer- 
tain. In fact, we do not even know whether an accretion disk 
forms at all. Agol & Kamionkowski (2002) show that most 
of BHs accreting gas should form accretion disks; but these 
disks are not necessarily thin. If the accreting gas is able 
to form a thin accretion disk (Shakura & Sunyaev 1974), 
then rj = 0.1, as assumed by Mii & Totani (2005). However, 
it seems to be more realistic that the gas forms an ADAF 
(i.e. Advection-Dominated Accretion Flow) disk, whose ra- 
diative efficiency^ is of the order of r; = 0.001 for IMBHs 
of mass M. ~ 10" M0 (Quataert & Narayan 1999). Finally, 



3 More precisely, equation (11) refers to the bolometric luminos- 
ity due to the Bondi-Hoyle accretion. However, detailed models of 
spectra of black holes accreting in the Bondi-Hoyle regime (Beskin 
&; Karpov 2005) or forming ADAF disks (Narayan, Mahadevan 
&; Quataert 1998) show that more than 60% of the bolometric 
luminosity is emitted in the X-ray range and more than 40% be- 
tween 0.2 and 10 keV (approximately the bandpass of Chandra 
and XMM). Because of the other uncertainties in our calculations, 
we think that the approximation that nearly all the Bondi-Hoyle 
luminosity is emitted in the X-ray band is acceptable. 



* We consider only the particular case of the equation by Mii & 
Totani (2005) in which M. > 1O3M0 

5 The luminosity for an ADAF disk scales as Af^, where M is 
the accretion rate. However, if log{M / M^dd) ~ — 4, —2, where 
M^dd is the Eddington accretion rate, the efficiency of an ADAF 
disk is about two orders of magnitude lower than the efficiency of 
a thin disk (see Figure 7 of Narayan ct al. 1998). We can assume 
that the efficiency of an ADAF disk is 77 = 10—3, because the 
accretion rates of the IMBHs we are considering fall in the range 
above. 
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models which take into account gas magnetization (Beskin 
& Karpov 2005) show that a high efhciency {rj « 0.1) is al- 
lowed, even if the thin disk does not form. To decide among 
these different models is beyond the scope of this paper; we 
will consider the two different values 77 = 0.1 and rj = 0.001 
bracketing the above possibilities in all our cases. 

In equation (12) fdisk is the fraction of IMBHs pass- 
ing through the molecular disk of the Galaxy, for which we 
assume a scale height zmc = 75 pc and a radial extension 
Rmc ~ 20 kpc (Sanders, Solomon & Scovillc 1984). Mii & 
Totani estimated fdisk 4.5 x 10"'', based on the hypoth- 
esis that IMBHs are a halo population following a standard 
NFW profile. Our simulations allow a more precise and di- 
rect estimate of fdisk from our simulations. As an example, 
in Fig. 2 we show a snapshot of our simulations, where the 
positions of IMBHs passing through the molecular disk are 
shown. Table 3 reports the simulated values for fdisk and 
NuLX with luminosity Lx > 10^® erg . 

For r/ = 0.1, we find that, if IMBHs arc born in 3 a fluc- 
tuations (corresponding to $7, « 10~^ni,) and their present 
distribution in the Milky Way follows the DMM model (case 
Al), the number of ULXs associated with these IMBHs is 
still consistent with zero [Nulx ~ 0.2; Table 3, third col- 
umn). Instead, if n. « lO^^ili, (case Bl), the Milky Way 
should host about 30 active ULXs. Then, we conclude that, 
if the IMBHs follow a DMM distribution, fi, « 10"^^^ can 
be considered as an upper limit for the present density of 
IMBHs. If, on the contrary, the IMBHs follow a standard 
NFW profile, as assumed by Mii & Totani (2005), the num- 
ber of ULXs obtained for f2, w 10~^f2i, (case B2) is still 
marginally consistent with zero. Finally, if IMBHs follow 
the DMM distribution but form only in fluctuations with 
a > 3.5, they arc so rare that no ULX is expected to be 
seen in the Milky Way (case CI). 

However, if IMBHs are surrounded by low efficiency 
ADAF disks (r) = 10"^), the upper limit for a DMM pro- 
file becomes O, = 0.1 ili,; whereas there are nearly no con- 
straints for the NFW profile. It is worth noting that Mii & 
Totani (2005) mainly considered the case of maximal effi- 
ciency (77 = 0.1), which is probably unlikely according to 
Agol & Kamionkowski (2002). 

Equation 12 tells us even another information: the dark 
matter halo of the Milky Way cannot be entirely composed 
by BHs with mass > 10=^ Mq. In fact, if M. = 10^ Mq and 
iV, = 10^ (corresponding to the assumption that the Milky 
Way dark matter halo is composed by BHs as massive as 
10^ M0), NuLx{> lO^^ergs-i) ~ 2 x 10^ r; for a DMM 
profile and Nulx{> lO^^erg s'^) ~ 9 x 10^ 77 for a NFW 
model. This means Nulx ^ 1 both for a thin and an ADAF 
disk model (unless rj -C 10~^). Then, the dark matter halo 
can be entirely composed by BHs only if their mass is less 
than 10* Mq , ruling out the scenario proposed by Lacey 
& Ostriker (1985), in which halo BHs can account for the 
galactic disk heating. 

3.2 Radial and luminosity distribution of ULXs 

The previous results can be refined by using our simula- 
tions instead of equation (12). In fact, from the simula- 
tions we know the number N{z < zmc) of IMBH particles 
which at a given time are in the molecular disk (defined 
by the scale height zmc = 75 pc and the radial extension 



Rmc ~ 20 kpc). It is well known that H2 in the Galaxy 
is not distributed uniformly within such disk, but it has a 
clumpy structure made of molecular clouds. Following Agol 
& Kamionkowski (2002) we derive the actual volume frac- 
tion of the molecular disk occupied by the clouds, i.e. 



/mc = 



(/?-2)(Emc) 



{P — 1)2 /imp ZMC rir, 
0.017, 



Tlmin / 



(1-/3)' 



(13) 



where (3 = 2.8 for an H2 cloud, (Emc) = 29 Mq pc'^ 
(Sanders et al. 1984; Mii & Totani 2005) is the average sur- 
face density of molecular clouds, mp is the proton mass, 
nmin = 10^ cm~^ and rimax = 10* cm~^ are the minimum 
and maximum density, respectively, observed in molecular 
clouds. Thus, the number of IMBHs which at a given time t 
arc embedded into a molecular cloud is N{z < zmc) f mo- In 
practice, wo randomly select from our simulations a fraction 
/mc of the IMBHs which at a given time t have z < zmc 
and R < Rmc- For this sample of IMBHs, we derive the 
Bondi-Hoyle luminosity Lx as in equation (11), assuming 
that Cs ~ 0.3 km s^'^ and ctmc ~ 3.7 km s^^ (Larson 1981; 
Solomon et al. 1987; Mii & Totani 2005). We then identify 
as ULXs those IMBHs which have Lx > 10^^ erg s"^ Aver- 
aging this number over the entire simulation, we obtain an 
estimate Nulx(> Lx) of the mimbcr of ULXs in the Milky 
Way. For all the considered cases, the rmmber Nulx(> Lx) 
(Table 3, fourth column), derived in this way, is consistent 
with the value Nulx{> Lx) (Table 3, third column), de- 
rived from equation (12), confirming the validity of the Mii 
& Totani calculation. 

This alternative method to derive the number of ULXs 
contains additional important pieces of information concern- 
ing the spatial distribution of ULXs and their luminosities 
(see Fig. 3 for the case Bl). These distributions are meaning- 
less for the Milky Way, where no ULXs have been detected. 
Nevertheless, it could be interesting to compare them with 
the distributions of ULXs observed in other spiral galaxies. 
Fig. 3 shows that, if a DMM profile is adopted for IMBHs, 
ULXs appear to be mostly concentrated towards the Galac- 
tic center. This seems to be at odds with observations, which 
have shown that ULXs of external galaxies tend to be pref- 
erentially located in spiral arms (Liu & Bregman 2005). We 
have to keep in mind, though, that the present calculation is 
based on the molecular hydrogen distribution of the Milky 
Way, which could be quite different from that of other galax- 
ies hosting ULXs; the latter are often starbursting, very gas 
rich systems. 

The predicted ULX luminosities (Fig. 3) are mostly in 
the range 1-5 xlO^® erg with only few sources show- 
ing luminosities higher than 10*" erg s~^. This rapid falloff 
of the number of ULXs for increasing X-ray luminosities is 
consistent with observations (Grimm et al. 2003). On the 
contrary, simulations following the accretion of IMBHs in 
binary systems indicate a number of low luminosity ULXs 
which is only a factor « 2 higher than the number of sources 
with Lx > 10"" erg s"^ (Madhusudlian et al. 2005). As 
a caveat, we recall that our calculations assume a Bondi- 
Hoyle luminosity, which might be a relatively oversimplified 
approximation. 
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Figure 3. Distribution of ULXs as a function of tlieir Galactocentric distance (left panel) and X-ray luminosity (right panel) for the 
case Bl and rj = 0.1, after O.bGyr. 



3.3 Non ULX sources 

By using the technique described iu Section 3.2, we can also 
derive the number of low luminosity X-ray (in brief, non- 
ULX) sources with Lx < 10^^ erg s"^ (Table 3; fifth col- 
umn). An interesting result is that, if rj = 0.1, IMBH lu- 
minosities are always as high as 10"^^ erg s^^ (sec Fig. 4, 
where the dotted line represents IMBHs accreting molecu- 
lar gas, including ULXs). Instead, if only ADAF disks can 
form, the luminosities reached by accreting IMBHs in molec- 
ular clouds are lower, in the range from lO'^'"" to lO''* erg 
s~^. Luminosities from 10'^'^ to 10"^^ erg s~^ are reached, in 
our Galaxy, only by a few supernova remnants (Vink 2006) 
and by high mass and low mass X-ray binaries (HMXBs, 
LMXBs; Psaltis 2006). Then, IMBHs accreting gas in molec- 
ular clouds should bo among the most powerful Galactic 
X-ray sources and therefore should have been already de- 
tected, provided they are not transient. A strong constraint 
on the density of IMBHs thus descends from the require- 
ment that the number of IMBHs with Lx < W'^^ erg s^^ is 
lower than the number of detected Galactic sources emitting 
at the same luminosities which lack of certain identifications 
with other kind of objects (such as HMXBs and LMXBs). 
This analysis will be carried out in Section 5, considering 
X-ray sources produced by IMBHs accreting both within 
molecular clouds and atomic hydrogen regions. 



4 IMBHS ACCRETING ATOMIC HYDROGEN 

Mil & Totani (2005) neglected in their analysis IMBHs pass- 
ing through atomic hydrogen regions, because their lower 
density ( ^ 1 cm~^) powers much lower X-ray luminosities 
than in molecular gas. (King ct al. 2001). However, atomic 
hydrogen is much more diffuse in the Milky Way than H2, 
and IMBHs are so massive that they can have non-negligible 
luminosity even accreting in such rarefied environment. Cur- 



rent models of the hydrogen distribution in the Milky Way 
(McKee & Ostriker 1977; Rosen & Bregman 1995) suggest 
the existence of three different phases: a neutral cold compo- 
nent (T ^ 10^ K), a warm (T « 10* K) and a hot (T ^ 10® 
K) component. In our work we neglect the hot component, 
since, even if its filling factor is high (up to 0.7, Rosen & 
Bregman 1995), it has am average density of w 10~^ cm~^ 
and a sound speed of « 100 km s~^, so that the X-ray lu- 
minosity of IMBHs accreting hot gas is expected to be very 
low. We define an atomic hydrogen disk as a disk having 
cut-off length Rh = 20 kpc (the data show an exponential 
fall of neutral hydrogen density outside 20 kpc; Lockman 
2002) and scale height zh = 100 pc (Baker & Burton 1975; 
Sanders et al. 1984; Dickey & Lockman 1990). Adopting the 
procedure described in Agol and Kamionkowski (2002), we 
derive the volume fraction occupied by cold neutral hydro- 



gen, fcH, as: 
fcH = 



{(3cH - 2)(Sch) 



(/3cfl" — 1)2 ^ nip ZH nmin,CH 
(l-/3cff)" 

_^ / rj /-> u \ 



,CH 



i,CH 



(14) 



where (3ch ~ 3.8 (Agol & Kamionkowski 2002), (Ech) is 
the average surface density of neutral hydrogen {{Ech) = 
4.5 Mq pc"^ if i? > 4 kpc and (Ech) w if i? < 4 kpc; Agol 
&: Kamionkowsky 2002; Sanders et al. 1984), rimin.CH and 
nmax,CH arc the minimum and maximum density of neutral 
hydrogen, respectively {rimin.cH ~ 1 cm"'^ rimax.CH ~ 5 
cm~^; Bregman, Kelson & Ashe 1993). Substituting these 
values into equation (14), we obtain fcH = 0.48 ii R > 
4 kpc and fcH « if J? < 4 kpc. This value is in good 
agreement with run E of Rosen & Bregman (1995), which 
is a suitable fit of cold and warm hydrogen observations 
(Dickey & Lockman 1990). For consistency, we assume that 
the filling factor of the warm component is fwH=0.2, as 
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Table 3. Results. 



Run 


J disk 


NuLx'' 


NuLx" 






Af(1036-39 erg s-^)f 


Al 


0.022±0.003 


0.2 (0.002) 


0.2±0.2 (0.002±0.002) 


1.0±0.6 (1.2±0.5) 


45±12 (45±12) 


18±7 (1.2±1.0) 


A2 


0.0008±0.0006 


0.008 (8x10-5) 


0±0 (0±0) 


0±0 (0±0) 


5±4 (5±4) 


0.4±0.4 (0±0) 


Bl 


0.025±0.002 


28 (0.28) 


40±6 (0.5±0.5) 


310±24 (350±24) 


4056±125 (4058±125) 


1650±70 (236±15) 


B2 


0.00090±0.00007 


1 (0.01) 


0.5±0.5 (0.007±0.007) 


5.6±1.5 (6.1±1.5) 


495±39 (495±39) 


148±21 (5±3) 


C 


0.039±0.018 


0.02 (0.0002) 


0±0 (0±0) 


0±0 (0±0) 


3±3 (3±3) 


0.09±0.09 (0±0) 



The values refer to a thin disk with rj = 0.1 (the values in parenthesis refer to an ADAF disk with 77 = 0.001). 
"Average fraction of IMBHs passing through the molecular disk (see Section 3.1). 

'Number of ULXs with Lx > 10'^^ erg s~^ derived from equation (12) adopting fj, = 2 (see Section 3.1). 
'^Number of ULXs with Lx > lO^^erg s~^ derived from our simulations (see Section 3.2 and Fig. 3). 

''Number of sources which accrete molecular hydrogen (see Section 3.3) with X-ray luminosities Lx < 10''^ erg s~^, derived from our 
simulations (see Fig. 4 and 5). 

"^Number of sources which accrete atomic hydrogen (see Section 4), derived from our simulations. All of them have Lx < lO'^^ erg s"^ 
(see Fig. 4 and 5). 

■^Number of sources which accrete atomic or molecular hydrogen and have X-ray luminosity 10'''' < Lx < lO'^^ erg s~^ (see Section 5). 





34 36 38 
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Figure 4. Distribution of all the X-ray sources as a function of the radial distance (left panel) and luminosity (right panel) for the case 
Bl and r] = 0.1, after 2 timescales. Open histogram: IMBHs passing through cold neutral hydrogen; light shaded histogram: IMBHs 
passing through warm hydrogen; heavy shaded histogram: IMBHs passing through molecular hydrogen. Although the distributions 
slightly change with time due to the small statistics, their main features remain unaltered. 



in run E of Rosen & Bregman (1995). Then, the number 
of IMBHs which, at a given time t, are passing through 
cold (warm) hydrogen regions is N{z < zh) fcH {^{z < 
zh) fwH)- Adopting the same procedure as in Section 3.2, 
we randomly select a fraction fcH (for cold hydrogen) or 
fwH (for warm hydrogen) of the IMBHs which, at a given 
time, pass through the neutral hydrogen disk, and we derive 
the Bondi-Hoyle luminosity® for each of them using equation 



^ We assume a turbulent velocity cth = 10 km s~^ both for cold 
and warm hydrogen (Lockman & Gehman 1991). The adopted 
sound speed is 1 km s~^ for cold hydrogen and 10 km s~^ for 
warm hydrogen regions. 



(11). The results are shown in Table 3, sixth column, and in 
Figure 4. 

If »7 = 0.1, IMBHs passing through cold neutral hydro- 
gen regions show luminosities of the order of IQ^'^"''^ erg s^^^, 
with a high luminosity tail at > 10^^ erg s~^; IMBHs passing 
through warm hydrogen regions produce lower luminosities, 
ranking from 10^^ to 10^^ erg s~^ with an extended high 
luminosity tail extending to 10"^* erg s~^ (Fig. 4). Due to 
the large filling factor of the atomic hydrogen with respect 
to the molecular one, the number of IMBHs accreting HI is 
a factor « 10 — 30 higher than the number of IMBHs which 
accrete H2, even if the luminosities are lower and nearly no 
ULX can be produced in atomic regions. 
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Figure 5. Distribution of all the X-ray sources as a function of their Galactocentric distance (left panel) and luminosity (right panel) 
for the case Bl and rj = 0.001, after 0.5 Gyr. Open histogram: IMBHs passing through cold neutral hydrogen; light shaded histogram: 
IMBHs passing through warm hydrogen; heavy shaded histogram: IMBHs passing through molecular hydrogen. Although the distributions 
slightly change with time due to the small statistics, their main features remain unaltered. 



Most interestingly, we note from the column 6 of Table 
3 that for Q, = 0.1 (both for r] = 0.1 and for rj = 0.001), 
the expected number of X-ray sources is huge, both for in 
the case of a DMM profile (case Bl, ~ 4000 sources) and for 
a NFW profile (case B2, « 500 sources). The reason the case 
B2 (which yielded the acceptable number of ULXs ~ 1) pre- 
dicts so many X-ray sources depends on the HI distribution 
properties: atomic gas is less concentrated than molecular 
clouds. Because a number of X-ray sources (not identified 
with HMXBs or LMXBs) > 500 is definitely too high for 
the Milky Way, we can robustly exclude Q, ^ 0.1 r^i,. 

For 77 = 0.001, such as for an ADAF disk, IMBH lumi- 
nosities are much lower, spanning 10^^"^^ to lO^'' erg s 
(Fig. 5). Even if the total number of sources remains nearly 
unmodified (Table 3; sixth column), the fact that most of 
them present luminosities <^ 10'^^ erg s~^ makes compari- 
son with observations more difficult. In the next section we 
attempt to constrain the Galactic IMBH density by compar- 
ison with X-ray observations, considering X-ray sources pro- 
duced collectively by IMBHs accreting both within molecu- 
lar clouds and atomic regions. 



5 COMPARISON WITH OBSERVATIONS 

In the previous Sections (3.1-3.2) we tried to constrain the 
density of IMBHs in our Galaxy by the fact that no ULXs 
have been detected in the Milky Way. From our simulations 
we found that some accreting IMBHs might also emit as 
non-ultraluminous X-ray sources, being in some cases bright 
enough to had been refiably detected by current X-ray satel- 
lites. 

Hereafter, we match the results of our simulations with 
the X-ray observations. In particular, we compare the pre- 
dicted IMBH X-ray emission with our knowledge of the X- 



ray sky in order to define an upper limit on the presence of 
these objects in our Galaxy. 

So far we have predicted IMBH X-ray luminosities (for 
a summary see Table 3) of ~ 10''^"''^ ergs~^, mostly de- 
pending on the assumed accretion efficiency rj, disk model 
and molecular or atomic accreting material. Searching in the 
observations for an upper limit of possible IMBHs in such a 
wide luminosity range is a non-sense, mainly because at the 
low luminosities many sources were certainly missed. 

What we then study here are only the IMBHs with a 
predicted hard X-ray luminosity between 10^® — 10^^ ergs"'^ 
(see Table 3, last column). In all these cases, the high lumi- 
nosity of these sources makes us confident that we should 
have seen them in the monitoring campaign of the Galaxy 
with the new generation satellites (within a certain distance 
depending of the flux resolution of the given satellite) . 

The most uncertain point is whether IMBHs accreting 
gas are transient or persistent sources. If the IMBH would 
be able to form a (thin or ADAF) accretion disk, it should 
also be a transient source. Instead, IMBHs accreting in the 
Bondi-Hoyle regime without forming a disk, as suggested by 
Beskin & Karpov (2005), should show flares; it remains un- 
clear if they can be transient sources or not. On the other 
hand, the IMBH could be transient also as a consequence of 
properties of the interstellar medium. In fact the accretion 
rate is roughly proportional to the density of the gas, and 
the scintillation measures show that the density fluctuations 
of the interstellar medium can be as high as a factor 100 on 
scales from ~ 10^** down to « 10^^ cm (Rickett 1990; Lam- 
bert & Rickett 2000; Cordes & Lazio 2001; Ferrara & Perna 
2001). A halo IMBH can easily travel « 10^^ cm in about 
one day, and thus could suffer, in principle, changes of a fac- 
tor ~ 100 in its flux in this range of time. As a consequence. 
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we have considered all the sources meeting our requirements, 
both persistent or transient during the observations. 

The most wide and sensitive survey available for our 
aims is the soft gamma-ray survey recently obtained by 
the IBIS/ISGRI (Imager on Board INTEGRAL Satel- 
lite/INTEGRAL Soft Gamma-Ray Imager) instrument on 
board of the INTEGRAL satellite (Bird et al. 2004, 2006). 
This survey observed 50% of the sky with a flux limit of 
1 mCrab in the 20-100 keV energy range. 

Among more than 200 sources detected by the 
IBIS/ISGRI soft gamma-ray survey scan, we excluded all 
the sources that certainly could not belong to the sample 
of possible IMBHs. In particular, wc excluded all the well 
established X-ray binaries, which are a well known highly 
luminous Galactic class (both as transient and persistent 
sources) . Furthermore we withdraw from our sample all the 
X-ray binaries known to host a neutron star (e.g. either be- 
cause showing pulsations or thermonuclear bursts) . Wo then 
filtered for a couple of highly energetic supernova remnants. 
After this first filtering we end up with a few tens of un- 
known objects. 

Given the fact that what IBIS/ISGRI measures is a cer- 
tain flux at Earth and not a luminosity, which is usually hard 
to derive because of the poorly known distances, we put the 
sample of sources we derived after the latter filtering, at dis- 
tances between 1-15 kpc , and we took all the sources with 
an inferred luminosity 10^® — 10^^ erg s^^, which implies in 
terms of flux all the uncatalogucd IBIS/ISGRI objects with 
a detected flux (within their errors) >4.8 mCrab in the 20- 
40 keV energy range. Note that this flux limit is derived plac- 
ing a source emitting 10'*' crgs^^ at 15 kpc (e.g. the edge of 
our Galaxy), it would be detected by IBIS/ISGRI at a flux of 
4.8 mCrab in the 20-40 keV energy range, well above the flux 
limit of the survey. Hence we are confident that, if present, 
our putative Galactic IMBHs would had been detected in 
the 50% of the Galaxy covered by the IBIS/ISGRI survey. 
Note that the flux limit of 4.8 mCrab we assumed includes, 
for the completeness of our analysis, the worst case of the 
faintest source at the largest distance: the fact that we are 
looking for an upper limit on the number of these possible 
IMBHs allow us to make this assumption. 

Under these assumptions, we found only 3 IBIS/ISGRI 
unidentified sources which match our requirements. These 
sources were all persistent during the IBIS observations. 

Their luminosity falls in the lO'^^-lO'''^ erg s~^ range, all of 
them close to the 10^*^ erg s"^ bound. As the IBIS/ISGRI 
catalogue covers 50% of the Galaxy, we then tentatively pre- 
dict an upper limit of 6 sources with these characteristics in 
the entire Galaxy, if the volume observed is a fair sample. In 
a few years all the Galaxy will be covered by the IBIS/ISGRI 
survey and our tentative extrapolation may be refined. 

Let us now compare this number with that of IMBHs 
predicted by our simulations in the same luminosity range 
and reported in the last column of Table 3. 

a) Thin disks 

For a thin disk, even case Al (fi, = 0.001^6, DMM 
profile) yields « 18 sources, a value three times higher than 
observed. Furthermore, 4 of these simulated sources have 
Lx > W^^ erg s~^. From an additional run with Q,, = 
10~*f2i, and the DMM profile (not reported in Table 3 for 
simplicity) we saw that the number of sources with 10^*^ < 



Lx < 10"*^ erg s^^ is 0.6±0.6. Wc conclude that the upper 
limit in the case of a Shakura-Sunyaev disk and a DMM 
profile is f2, = 10~* — 10~^f2i,, similar to the upper limit 
found by the number of ULXs alone (see Section 3.1-3.2). 
Instead, for a NFW profile the allowed density of IMBHs is 
> f2, = 10~^f2i, (case A2; corresponding to 0.4±0.4 sources), 
but definitely < f2, = 10"^ fit (case B2; 148±21 expected 
sources), strengthening the constraint we found from the 
number of ULX. 

b) ADAF disks 
For the more realistic case of an ADAF disk, the 
constraints we obtain from the comparison with the 
IBIS/ISGRI sources arc stronger than for the number of 
ULXs alone. In fact, if wc assume a DDM profile, the up- 
per limit for the density of IMBHs is about n. = 10~''O(, 
(case Al; 1.2±1.0 expected sources), much lower than f2, = 
10~^f26 (case Bl; 236±15 expected sources), derived from 
the number of ULXs. If we consider a NFW model, the up- 
per limit is f2, = 10~^Oj, (case B2; 5±3 expected sources); 
whereas there were no significant constraints from the ULXs. 
In summary, we must take with care the results of this com- 
parison between simulated and observed X-ray sources, be- 
cause of the huge uncertainties of our model. However, from 
the comparison with the IBIS/ISGRI unidentified sources 
we derive, in general, much stronger constraints than from 
the number of ULXs. 



6 CONCLUSIONS 

In this paper we have simulated the dynamical and emis- 
sion properties of putative IMBHs which could inhabit our 
Galaxy. IMBHs arc modeled as a halo population, dis- 
tributed following a NFW or a more concentrated DMM 
profile. We assumed that IMBHs, passing through molec- 
ular or atomic hydrogen regions, could accrete gas, form- 
ing X-ray sources (either ultra- luminous or not). From the 
comparison of our simulations with the number of ULXs in 
the Galaxy (Section 3.1-3.2) and with the non-ultraluminous 
unidentified X-ray sources in the IBIS/ISGRI catalogue, we 
have derived the most stringent (to our knowledge) upper 
limits on the density Q, of IMBHs. The main results can be 
summarized as follows: 

• If IMBHs accrete with efficiency rj = 0.001 (i.e. via an 
ADAF disk), we obtain a strong upper limit f2, < 10~^n6 
for a DMM profile and Q, = IQ-^Qb for a NFW profile. 

• If the IMBHs accrete with efficiency rj — 0.1 (i.e. if a 
thin accretion disk around the IMBH is formed), the upper 
limit of the density of IMBHs is D,. = lO'^fit for a DMM 
profile and Q. = lO'^nj, for a NFW profile. 

These results are still affected by some model uncertain- 
ties, as the emission mechanism and the IMBH distribution. 
In addition, computational requirements have forced us to 
use high and equal mass (m. = 5 x lO* Mq) IMBH parti- 
cles. Constraints for lower IMBHs masses are expected to 
be weaker. We can guess how the above upper limits change 
for different values of the IMBH mass by using the equa- 
tion (12). For example, if wc assume rn, = 77 = 0.1 
and a DMM profile, the upper limit of the IMBH density 
becomes f2, = 10~^Qb, about one order of magnitude lower 
than for m, = 10* Mq. As a further caveat, this extension to 
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lower masses is possible only for the comparison with ULXs 
(and not with IBIS/ISGRI sources), because it is based on 
eq. (12). Therefore, higher resolution simulations would be 
required to extend our studies to lower mass IMBHs or to 
consider a more realistic IMBH mass spectrum. Higher res- 
olution simulations (where the mass of star particles can be 
orders of magnitude lower than the mass of IMBH particles) 
are also needed to take into account the dynamical friction, 
which could play a crucial role. 

Another caveat concerns the validity of the DMM pro- 
file. The simulations by DMM neglect the contribution of 
IMBHs in building up SMBHs, either by mergers (Islam, 
Taylor, & Silk 2003, 2004a,b,c) or by accretion and close 
dynamical encounters (Volonteri et al. 2003). Monte Carlo 
simulations combined with semi-analytical models (Volon- 
teri & Perna 2005) show that, if all these factors axe taJsen 
into account, the number of IMBHs could be up to 2 or- 
ders of magnitude lower, leading to an f2, ~ 10 — 100 lower 
than our estimates, and therefore compatible with the non- 
detection of ULXs in the Milky Way. However DMM take 
into account the bias in the formation sites of IMBHs, the ac- 
cretion into larger halos, the role of both dynamical friction 
and tidal stripping, which were neglected or described by 
rough models in the previous studies. Unfortunately current 
simulations cannot account for all these effects at the same 
time. In conclusions, even if our results could be improved 
under many aspects, we consider it as a success that our 
models strengthen by a factor 10-1000 the currently adopted 
upper limits for the density of IMBHs (i.e. O, « 0.02; van 
der Marel 2004). 
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